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ABSTRACT 

We observed PSR B1259-63, a young non-accreting pulsar orbiting around a Be star SS 2883, eight times 
with the Suzaku satellite from July to September 2007, to characterize the X-ray emission arising from the 
interaction between a pulsar relativistic wind and Be star outflows. The X-ray spectra showed a featureless 
continuum in 0.6-10 keV, modeled by a power law with a wide range of photon index 1 .3-1.8. When combined 
with the Suzaku PIN detector which allowed spectral analysis in the hard 15-50 keV band. X-ray spectra do 
show a break at ^ 5 keV in a certain epoch. Regarding the PSR B 1259-63 system as a compactified pulsar 
wind nebula, in which pairs are assumed to be accelerated at the inner shock front of the pulsar wind, 
we attribute the X-ray spectral break to the low-energy cutoff of the synchrotron radiation associated with 
the Lorentz factor of the relativistic pulsar wind 71 ~ 4 x lO*". Our result indicates that Comptonization of 
stellar photons by the unshocked pulsar wind will be accessible (or tightly constrained) by observations with 
the Fermi Gamma-ray Space Telescope during the next periastron passage. The PSR B 1259-63 system allows 
us to probe the fundamental properties of the pulsar wind by a direct means, being complementary to the study 
of large-scale pulsar wind nebulae. 

Subject headings: acceleration of particles — radiation mechanisms: non-thermal — pulsars: individual (PSR 
B1259-63) — X-rays: binaiies 



1. INTRODUCTION 

PSR B 1259-63 is a young radio pulsar (spin period 48 
ms) orbiting a B2e star SS 2883 in a highly eccentric 3.4 

i T orbit, having spindown power of iSp :^ 8 x 10 ^^ erg s"' 
Johnston et all 1992i[T994t [Manchester et al.ll995h. The dis- 
tance from the Earth is estimated as ^ 1 .5 kpc (iJohnston et aTl 
11994 . PSR B 1259-63 was the first example of a radio 
pulsar forming a binary with a non-degenerate companion 
( [Johnston et"aLl il992). The PSR B1259-63 system is an im- 
portant astrophysical laboratory for the study of the pulsar 
wind interacting with a stellar wind, and the nonthermal ra- 
diation arising from the shocked pulsar wind. Apart from 
PSR B 1259-63, the only other example k nown in our Galaxy 
is the PSR J1740-3052 binary system (IStairs et all 1200 ih . 
which has, however, much lower spindown power 

Be stars are characterized by strong mass outflows 
through the formation of a slow and dense equatorial disk 
dWaters et al. 1988). Because PSR B1259-63 has a highly 
eccentric orbit and the equatorial disk is likely inclined with 
respect to the orbital plane, the pulsar is expected to inter- 
act with the disk of the Be star near periastron, by crossing 
the disk tw ice before and after the periastron passage (e.g., 
iMelatos eT al. 1995). The signatures of disk crossings can be 
seen in the radio observations, such as pulsar eclipse, the light 
curve of transien t unpulsed emission, and th e change of dis- 
persion measure (IJohnston et al.lfT99a.l2005h . The disk cross- 
ings also play important roles in producing the X-ray and 
gamma-ray emission (Chernyakova et al. .2006i) . 
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Enhanced high-energy emission around periastron which 
is thought to arise from the interactions between the 
pulsar wind and the stellar win d /disk, has been ob- 
served in X-ray (KasDi et al.' 119951: iHirayama e t aT 19961 
1999; Chernvakova et al. 200^), soft 7-ray ([Grove et al] 
[l995; Shaw et al. 2004), and very high energy gamma-rays 
dAharonian et al.. .20051) . The basic properties of the X-ray 
emiss ion have been revealed by the series of ASCA observa- 
tions ("Kasp i et al.|[l995t IHirayama etaH 1 19961 [1999), which 
are summarized as follows: (1) The X-ray luminosity is mod- 
est around periastron with ^ 10^'^ erg s"', corresponding 
to about one percent of the spindown power. The intensity 
within a day is constant, while it varies by a factor of ~ 2 dur- 
ing the periastron passage. At around apastron, the luminosity 
decreased by one order of magnitude; (2) The X-ray spectrum 
is characterized by a power law of photon index F ~ 1 .6-2.0 
without any detectable line emission; (3) No significant X- 
ray pulsed emission at the pulsar spin period has been found. 
The upper limit (at a 90% level) on the pulse component is 
~ 10% of the total X-ray flux. Observations during the pe- 
riastron passage in 2004 with XMM-Ne wton have revealed 
more detailed behavior of X-ray emission (IChernyakova et afl 
[2OO6I) . The X -ray brightening prior to periastron was roughly 
coincident with the brightening of unpulsed radio emission, 
which can be identifie d with the entrance t o the Be star equa- 
torial disk. Moreover, IChernyakova et all ([2006) have found 
remarkable hardening down to F ~ 1 .2 at the entrance of the 
disk prior to periastron, followed by a recovery to F ~ 1 .5 
during the disk crossing. 

In the hard X-ray domain, a dedicated three-week CGRO 
observation in 1994 has resulted in a significant detection in 
the 30-200 keV range at a few mCrab level with the OSSE 
instrument, with the spectrum being cha racterized by a power 
law of F = 1 .8 ± 0.6 (Grove et al.lll995h . The OSSE spectrum 
was consistent with the time-averaged A5CA spectrum extrap- 
olated to the OSSE band. With the INTEGRAL satellite, the 
hard X-ray emission in the 20-200 keV band was detected 
with ^ 4(7 significance based on the observation performed in 
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21-25 March 2004 dShaw et alj|2004 . 

The properties of the X-ray and soft 7-ray emission ob- 
served with ASCA and CGRO have been shown to be con- 
sistent with the idea that the emission is produced by shock- 
accelerated pairs as a result of the interactions between the 
pulsar wind and the Be star outflows (iTavani. Arons. & Kaspil 
ll994tlKaspi et alJlT99llT^vani & Aronsll 19971) . The pairs 
are considered to be accelerated at the inner shock front of 
the pulsar wind and adiabatically expanding in the relativistic 
flow of the pulsar cavity. As a scale-down version of pulsar 
wind nebulae (PWNe), the PSR B 1259-63 system can be re- 
ferred to as a compactified pulsar wind nebula, and serves as 
a unique diagnostic of shock acceleration subject to very fast 
cooling. Indeed, synchrotron radiation by the accelerated 
pairs offers an excellent explanation for the observed X-ray 
emission, indicating the existence of a remarkably efficient 
acceleration mechan ism at an energy of Te V on a timescale 
less than ^ 100 sec ( iTavani & Arons|[l997h . The anisotropic 
inverse-Compton scattering on the intense stellar light of the 
accelerated pairs \ yas pred icted to be detectable by air 
Cherenkov telescopes dKirk et al. 1999). 

During the periastron passage in 2004, very high energy 
gamma -rays have been indeed detected by the HESS tele- 
scopes dAharonian et alj |2005). The level of the TeV gamma- 
ray flux and the spectral shape were in a reasonable agree- 
ment with the theoretical expectation based on the compact- 
ified PWN model. However, the light curve around perias- 
tron disagreed with the early expectation , requiring a revision 
of the model^. iKhangulyan et alJ d2007h have introduced the 
change of the maximum energy of pairs as a viable explana- 
tion for the behavior of the TeV light curve. Also, they argued 
that the increase of the adiabatic loss rate near periastron is 
another possible way to naturally account for the reduction of 
the TeV flux. Since the change of the maximum energy and/or 
the adiabatic loss rate can be imprinted in the synchrotron X- 
ray spectrum, the detailed X-ray observations are thus quite 
important to un derstand the behavior of the TeV radiation as 
emphasized bv Khangulvan et akl d2007h . 

In this paper, we present the results from our monitoring 
observations made with the Suzaku satellite in the X-ray and 
hard X-ray bands, with an effective energy range of 0.6-50 
keV, performed during the most recent periastron passage in 
2007. Combined with the previous measurements in X-ray, 
the Suzaku observations have clarified the spectral evolution 
of the PSR B 1259-63 system. Thanks to the wide band cov- 
erage, we have discovered the existence of a spectral break 
in a certain epoch. This paper is organized as follows. The 
observations and basic data reduction are described in ^ 
The results of X-ray analyses are presented in ^ In ^ 
we outline our model of a compactified pulsar wind nebula, 
and interpret the X-ray spectral break as a manifestation of 
the Lorentz factor of the pulsar wind. Throughout this paper, 
we adopt the distance to t he PSR B 1259-63 sy stem from the 
Earth to be D = 1.5 kpc (iJohnston et alJll994t) . For the cal- 
culation of the pulsar orbit, we assume a mass of SS 2883 to 
be Mc = IOMq and the pulsar mass to b e Mp = 1 AMp), whic h 
implies an inclination angle of / ~ 36° d Johnston et al.lfT994h . 
The adopted orbital elements of the pulsar are an orbital pe- 
riod of 1236.72 day, eccentricity of e = 0.8699, longitude of 
periastron of a; = 138.7°, and the epoch of periastron of MJD 

^ Even departing from the framework of the compactified PWN, a 
hadronic origin of TeV gamma-rays ( Kawachi et al. 2004) may replace the 
synchrotron-IC model (see also lChemvakova et alj2009) . 



48124.34 dManchester et ai][T995l) . AU errors are quoted at a 
90% confident level unless otherwise stated. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed the PSR B 1259-63 system eight times with 
the Suzaku satellite each with ~ 20 ks exposure in 2007 July, 
August, and September. The total exposure time amounts to 
~ 170 ks after standard data screening. We note that this is 
the first and possibly last chance for Suzaku to observe the 
first disk transit since target visibility does not permit us to 
observe it next time. Table[T]gives the log of the Suzaku obser- 
vations; the eight observing epochs are referred to as Szl-Sz8 
hereafter During the first four epochs (Szl-Sz4) simultane- 
ous observations with HESS were coordinated. The results of 
the HESS observations will be presented separately. 

The Suzaku observations were performed with the X-ray 
Imaging Spectrometer (XIS: iKoyama et alJl2007h in 0.3-12 
keV and the Hard X-ray Detector (HXD: iTakahashi et al] 
I2OO7I) in 13-600 keV. The XI S, located at the focal plane 
of the X-ray telescopes (XRT: Serlemits os et aljr2007h . con- 
sists of one back-illuminated CCD camera (XIS-1) and three 
front-illuminated CCDs (XlS-0, -2 and -3). One of the front- 
illuminated CCDs, XlS-2, was not available at the time of our 
observations, since it suffered from a fatal damage on 2006 
November 9, and unusable since then. The HXD consists of 
the silicon PIN photo diodes (hereafter PIN) capable of obser- 
vations in the 13-70 keV band and the GSO crystal scintilla- 
tors (hereafter GSO) which cover the 40-600 keV band. 

The XIS instruments were operated in a normal full-frame 
clocking mode, with a frame time of 8 s. In this work, we 
present results from front-illuminated CCDs (XlS-0 and -3), 
which have a larger effective area at high energies, and there- 
fore suited for our purpose. XlS-0 and -3 have almost iden- 
tical properties, so we co-added the data to improve photon 
statistics. PSR B1259-63 was observed slightly off-axis (3.5' 
from the optical axis), placed on the HXD nominal pointing 
position, to maximize the effective area of the HXD detectors. 
The HXD instruments were operated in the normal mode. The 
relative normalization of the HXD-PIN spectrum to the XIS 
was fixed at 1.15, as was determined by the calibration obser- 
vations of the Crab Nebula (Kokubun et al. 2007). 

The field-of-view (FOV) of the XIS contains another X-ray 
source, 2RXP J130159.6-635806 (shortly 2RXP J1301), as 
shown in Figure [T] This i s an accretion-powered p ulsar with 
a spin period of ^ 700 s (Chernyakova et af]l2005h . which is 
located 10' away from PSR B 1259-63. We have monitored 
the 700-sec accretion-powered pulsar with the XIS by tuning 
a telescope roll angle, to estimate its contribution to the hard 
X-ray flux measured by the non-imaging HXD-PIN detector 

We analyzed data delivered to us via pipeline processing 
version 2.0.6.13 (Szl-Sz5) or version 2.1.6.15 (Sz6-Sz8). 
Basic analysis was done using the HEASOFT software pack- 
age (version 6.5), with calibration files distributed on 9 July 
2008. For the XIS and HXD-PIN, we made use of cleaned 
event files, in which standard screening was applied. On 
the other hand, we reprocessed the HXD-GSO data by our- 
selves to apply appropriate gain files. The standard screen- 
ing procedures include event grade selections, and removal 
of time periods such as spacecraft passage of the South At- 
lantic Anomaly (SAA), intervals of small geomagnetic cutoff 
rigidity (COR), and those of a low elevation angle. Specif- 
ically, for the XIS, the elevation angle larger than 5° above 
the Earth and larger than 20° from the sunlit Earth limb are 
selected. The data screening criteria for the HXD include 
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PSR B1 259-63 



2RXPJ1301 59.6-635806 



Fig. 1. — Suzaku XISO 1-10 keV image (epoch Szl), shown with 
square root scaling. The FOV of the XIS is 17.8' x 17.8'. An 
accretion-powered pulsar 2RXP J130159.6-635806 is located 9.6' 
away from PSR B 1259-63. 



TABLE 1 
Log of Suzaku Observations 



ID 


Date 


MJD 


r 


Exposure (ks) 








(days) 


XIS 


PIN 


Szl 


2007 07 07 


54288.6 


-19.3 


21.9 


25.3 


Sz2 


2007 07 09 


54290.7 


-17.2 


19.5 


25.4 


Sz3 


2007 07 1 1 


54292.6 


-15.3 


22.7 


22.9 


Sz4 


2007 07 13 


54294.7 


-13.2 


22.9 


19.7 


Sz5 


2007 07 23 


54304.3 


-3.6 


19.7 


16.7 


Sz6 


2007 08 03 


54315.3 


4-7.4 


24.0 


20.1 


Sz7 


2007 08 18 


54330.1 


4-22.2 


20.5 


18.1 


Sz8 


2007 09 05 


54348.2 


-H40.3 


18.3 


20.8 



Note. — Modified Julian Date (MJD) and days 
from periastron (r) refer to the start time of the ob- 
servation. The exposure time is the net integration 
time after standard data screening. 



TABLE 2 

Results of PSR B 1259-63 Suzaku XIS Spectral Fitting 



ID 




r 


^"1-10 


xJ(d.o.f.) 




10^2 cm-2 


10- 


erg cm - s ' 




Szl 


0.50 ±0.02 


1.64 ±0.02 


25.8 ±0.3 


1.06 (354) 


Sz2 


0.51 ±0.03 


1.58 ±0.03 


26.7 ±0.3 


1.04(309) 


Sz3 


0.51 ±0.03 


1.35 ±0.03 


21.6±0.3 


1.19 (272) 


Sz4 


0.52 ±0.03 


1.44 ±0.03 


23.1 ±0.3 


1.30 (298) 


Sz5 


0.50 ±0.03 


1.83 ±0.03 


19.0±0.2 


0.85 (249) 


Sz6 


0.52 ±0.03 


1.73 ±0.03 


12.9 ±0.2 


0.92 (205) 


Sz7 


0.50 ±0.02 


1.69 ±0.02 


34.4 ±0.4 


1.01 (402) 


Sz8 


0.46 ±0.03 


1.57 ±0.03 


22.7 ±0.3 


1.12(253) 



Note. — Fitting Suzaku XIS spectrum of PSR B 1259-63 by 
a power law with photoelectric absorption in 0.6-10 keV. Ab- 
sorbing column density A'h, photon index F, and the 1-10 keV 
flux Fi_io (not corrected for absorption) are shown with 90% er- 



rors. 

the COR larger than 6 GV, and the elevation angle above the 
Earth limb greater than 5°. 

3. ANALYSIS AND RESULTS 

3.1. Power-law fit to XIS spectra 

The XIS spectra were extracted from circular regions of 3' 
radius centered on the source. Background spectra were ex- 



tracted from source free regions. Response matrices ("rmf") 
were generated by xisrmf gen, and ancillary response func- 
tions ("arf") were simulated by running xissimarf gen. 
The source count rate is typically 1 counts s"' per XIS in 0.6- 
10 keV, while the background count rate is only about 3 x 10"^ 
counts s"'. Source spectra were binned to a minimum of 100 
counts per bin. We performed spectral fitting in a range of 
0.6-10 keV after excluding a narrow bandpass around the Si- 
K edge (1.8-1.9 keV) where some calibration errors remain. 

Given that no previous observations have reported spec- 
tral changes within a day, we first performed model fitting of 
epoch-by-epoch X-ray spectra of the PSR B 1259-63 system 
accumulated for ~ 20 ks. As expe cted from past X-r ay ob- 
servations (e.g.. lKaspi et a l. 1995; Chernvakova et al.] |2006). 
all the XIS spectra of the PSR B 1259-63 system show a fea- 
tureless continuum without any emission lines. The Suzaku 
X-ray spectra are shown later in Figure |6l The XIS spec- 
trum (in terms of number flux per energy) was modeled by 
a power law, namely F(e) = Ke~^ where K is the normal- 
ization and r is the photon index, attenuated by an absorp- 
tion factor of exp[-A^Hc(e)]. Here A^h is the absorption col- 
umn density and a{e) represents the photoelectric absorption 
cross section. The results of the power-law fit are summa- 
rized in Table |2l The best-fit photon indices are found in a 
wide range of F = 1.35-1.83, while the absorption column 
density in a narrow range of A^h = (0.46-0.52) x 10^^ cm"^. In 
most cases we obtained a statistically acceptable fit. However, 
in the case of Sz4, the absorbed power-law fit was unaccept- 
able with a reduced chi-square of = 1-30 for 298 degree of 
freedom (d.o.f.). Also, the fitting quality for Sz3 seems bad 
with = 119 for 272 d.o.f. It should be noted that a quite 
hard power law was obtained in both epochs. Inspection of 
the fit residuals suggests the presence of spectral steepening. 
Thi s iss ue shall be investigated in detail using the HXD-PIN 
in El 

Figure |2] plots the flux measured by the Suzaku XIS as a 
function of t (days from periastron), together with a compila- 
tion of previous measurements. The new Suzaku data points 
largely help us track the time evolution of the X-ray flux from 
the PSR B 1259-63 system, and make it clear that the stel- 
lar disk plays essential roles in regulating the X-ray emis- 
sion from the PSR B 1259-63 system. In the X-ray lightcurve 
as a function of r (Fig. |2|i, there appear two bumps, which 
can be identified with the disk passages of the pulsar. The 
position of the post-periastron peak at epoch Sz7 (t = +22 
days) coincides with the enhancement of the TeV flux mea- 
sured with HESS in 2004, as well as the peak in the ra- 
dio lightcur ves. In accordance with the disk geometry in- 
troduced by IChemvakova et alj (l2006h . we associate this fea- 
ture with the pulsar's re-entrance to the Be star equatorial 
disk. The relatively smooth multi-orbital lightcurve over a 
decade suggests that orbit-by-orbit differences in X-ray be- 
havior would be small, though the current data sets are not 
enough to draw a definitive conclusion. We note that the un- 
pulsed radio lightcurves show significant differences from one 
periastron to next (see Fig. |2]i, which would support the idea 
that t he radio el ectrons originate from the shocked equatorial 
disk (.Ball et al.iri999) . 

In Figure |3] the best-fit values of F and Ah are shown 
as a function of r. As mentioned above, the column den- 
sity is found to be almost constant. Ah — 0.5 x 10^^ cm"^, 
over the course of Suzaku monitoring. This is roughly con- 
sistent with the values obtained with XMM-Newton during 
the previous disk passage in 2004, Ah — 0.45 x 10^^ cm 



Uchiyama et al. 



^ 3 — 



> 

<B 
O 

T 



2 — 



1 — 



> 

<D 



OO 



O 

A 



N 

o 



1.2§ 

1.00 

0.75 

0.50 

0.25 

0.00 
60 

40 

20 



ASCA 94 O 
SAX 97 □ 
XMM 04 O 
Suzaku 07 # 




. Radio 19197 (solid) 
Radio 2000 (dashed) 
Radio 2004 (dotted) 




-50 



50 
T (Days from Periastron) 



100 



Fig. 2.- 



X-ra- 
10"" erg cm"^ 



TeV, and radio fluxes as a function of time from periastron passage. (Top) 1-10 keV flux lightcurve in units of 
Absorption is not corrected for. Errors are s mall compared w ith the mark size. The Suzaku results obtained in this 
work are drawn as filled circ les. The ASCA GIS data points in 1994 ll Hiravama et alj[l996) . the Beppo-SAX ones in 1997, and the XMM ones 
in 2004 jChemvakova et an [2006) are shown as open triangles, squares and circles, res pectively. The periods ov er which the pulsar crossed 
the equatorial disk are shown (blue hatched regions) based on the geometry proposed by Chemvakova et al] ( l2006h . {Middle) Gamma- ray flux 
above 0.38 TeV in uri its of 10"" cm"^ s"' obtained with HESS jAharonian et alj2005h . {Bottom) 1.4 GHz radio lightcurves in 1997, 2000, 2004 
dJohnston et al.ll2005h . 



jChernyakova et all 120061) . It should be noted that smaller 
column density (A^h ~ 0.3 x 10^^ cm"^) has been observed 
with XMM-Newton when the pulsar was positioned outside 
the disk (Chernyakova et al. 2006). The excess column den- 
sity can therefore be ascribed to the dense stellar disk. On the 
other hand, photon index exhibits remarkable changes; sig- 
nificant spectral flattening characterized by F ~ 1.35 was ob- 
served during the first disk transit (r ^ -15 d). Such spectral 
hardening was discovered in 2004 with XMM-Newton, but it 
appeared slightly earlier in phase (r ^ -20 d). Later in sec- 
tion 33.41 we will demonstrate the presence of a spectral break 
during the hard-spectrum state, by utilizing hard X-ray mea- 
surements in the 15-60 keV band with the Suzaku HXD-PIN. 

3.2. Lightcurves 

The X-ray lightcurves for individual observations were ex- 
amined in order to study the temporal characteristic of the X- 
ray emission on short time scales. The 1-10 keV lightcurves 
with a bin size of 500 sec are shown in Figure |4] We found 
moderate hour-scale flux variability {^ 30%) at the time of 
disk crossings. On the other hand, the lightcurves in epochs 
Sz5 and Sz6 can be fit well with a constant, suggesting that 
the flux variations would be characteristics of disk passage. 
The size of the emitting regions may be smaller in the equa- 
torial disk. Alternatively, there is a large density fluctuation 
inside the disk. 

To search for possible spectral changes within each obser- 
vation, we subdivided each epoch into several periods typi- 
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Fig. 3. — Suzaku XIS best-fit parameters as a function of time 
from periastron passage. Statistical errors are smaller than the size of 
marks. {Top) The same lightcurve as the top panel of Fig.|2] {Middle) 
Photon index P. {Bottom) Absorption column density A^h- 

cally with an hour exposure, and fitted the short time-frame 
spectra with a power law modified by photoelectric absorp- 
tion. Though the power-law normalization differs between 
sub-epochs, no statistically significant spectral changes were 
found, which indicates the stability of a spectral shape on hour 
scales. 

The density fluctuation inside the Be star disk would be im- 
printed in the chaotic behavior of the X-ray lightcurves on 
hour scale in one way or another. The hour-scale variabil- 
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ity Suzaku uncovered would offer a way to infer the pulsar's 
location relative to the Be disk plane. Future intense cover- 
ages during the post-periastron disk transit of the pulsar with 
Suzaku will be able to constrain the structure of the Be star 
disk. 
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Fig. 4.— Suzaku lightcurves of PSR B1259-63 in the 1-10 keV 
band. The count rate with a bin size of 500 s from the start time 
of the observation is shown without background subtraction, whose 
contribution is only a few percent. 



3.3. Preparation for Broadband Analysis 

With the Suzaku observations, we are now accessible to a 
day-by-day broadband X-ray spectrum spanning two orders of 
magnitude in energy from 0.6 to 50 keV with sufficient sen- 
sitivities, by combining the XIS (0.6-10 keV) and the HXD- 
PIN (15-50 keV) instruments. 

To extract the HXD-PIN spectra of PSR B 1259-63, we 
have to subtract the following background spectra from the 
total cleaned spectra: (1) instrumental background so-called 
non X-ray background (NXB), (2) cosmic X-ray background 
(CXB), that is the sum of unresolved active galactic nuclei in 
the FOV of the PIN, and (3) the emission from the nearby 
pulsar 2RXP J1301 that Hes within the FOV of the PIN. 
The instrumental background spectra (perhaps due mostly 
to atmospheric neutrons), provided by the HXD instrument 
team, are constructed by the empirical models, which rely 
on a database of Earth occultation observations. The HXD 
team provides two kinds of background models, "quick" and 
"tuned". We found the two models give almost identical re- 
sults within the statistical uncertainties. In this paper we make 
use of the "tuned" NXB model, which offers better back- 



ground reproducibility ach^ the la systematic error of 
<3% dFukazawa et alj 120091) . The systematic error of the 
NXB subtraction does not have noticeable impact on the re- 
sults presented in this paper. 




5 10 
Energy [keV] 

Fig. 5. — Suzaku XIS and PIN spectra obtained by combining 
epochs Sz5 and Sz6. The XIS spectrum of PSR B 1259-63 is shown 
in blue, and that of 2RXP J 1301 in red. The PIN spectrum (black 
crosses) after background subtraction is decomposed into the two 
components through joint fitting: PSR B 1259-63 power-law model 
(blue line), and 2RXP 11301 cutoff power-law model (red line). 

The CXB flux in the FOV of the PIN instrument was 
subtracte d from each PIN spectrum in the same way as 
iTakahashi et al.i (l2008h . The CXB flux amounts to - 5% of 
the NXB, comparable to the systematic error of the NXB 
model itself. 

We now describe the procedures used for the derivation 
of hard X-ray spectra of 2RXP J1301 in each pointing. In 
fact, the main contributions to the systematic error of the 
PIN spectral data of the PSR B1259-63 system itself would 
be the uncertainties in the hard X-ray flux contaminated by 
2RXP J 1 30 1 . To describe the X-ray spectrum of 2RXP J 1 30 1 , 
we adopt the following empirical spectral form that has been 
traditionally applied to accretion-powered pulsars: 



F(e) = 



(1) 



Ke-' e < Ec 

Ke~' exp (-( e - Ec) / Cf ) e > Ec , 

where K is the normalization (IWhite et al.llI983h . The index s 
was frozen to be s = I , a typical v alue for accretion-powered 
X-ray pulsars (e.g.. lNagasellI989h . to reduce the number of 
free parameters. Noncyclotron puls ars can be characterize d 
by =4-15 keV and ef= 10-30 keV dMakishima et al.lI999l) . 

We chose the epochs Sz5 and Sz6 for the use of deter- 
mination of the spectral shape of 2RXP J1301 in the 1-50 
keV band, because the relative contribution of 2RXP J 1301 
in a PIN spectrum would be large in the two epochs, where 
the PSR B 1259-63 emission is relatively weak and steep. 
To obtain the model parameters relevant for 2RXP J 1301, 
the XIS spectrum of PSR B 1259-63 in 0.6-10 keV, that of 
2RXP J1301 in 1.5-10 keV, and the PIN spectrum in 15-50 
keV after subtraction of the background (NXBh-CXB) were 
jointly fit by a two component model consisting of a power 
law that describes PSR B 1259-63, and equation Q that de- 
scribes 2RXP J 1301. Each spectral component is attenu- 
ated by photoelectric absorption. The spectral data for the 
two epochs were summed to reduce statistical errors, because 
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the XIS spectra of PSR B 1259-63 (and also 2RXP J 1301) 
are similar in the two epochs. An energy-independent trans- 
mission factor (about 0.7) of the PIN collimator at the posi- 
tion of 2RXP J 1301 was taken into account when perform- 
ing a joint fit. The assumption of a power-law spectrum for 
PSR B 1259-63 in the 0.6-50 keV band in epochs Sz5 and 
Sz6 can be Justifi ed, to some extent, by the OSSE results 
jGrove et al.lll995l) . which came from the three-week obser- 
vation covering a certain orbital phase that encompasses the 
two Suzaku epochs. 

The joint fit was statistically acceptable with a reduced chi- 
square of xi = 1-02 for 768 d.o.f. In Figure|5]the XIS spectra 
of both PSR B 1259-63 and 2RXP J 1301, and the PIN spec- 
ti-um (the sum of PSR B 1259-63 and 2RXP J 1301 contribu- 
tions) are shown together with the best-fit model. We ob- 
tained the following parameters for the 2RXP J1301 model: 
Nh = (2.62 ± 0.08) X 10^2 cm'^, = 4.4!° j keV, and ef = 
14. 1 ± 1 .2 keV. The values seem reasonable for noncyclotron 
pulsars. The best-fit parameters for the PSR B1259-63 model 
are A^h = (0.50 ±0.02) x 10-^ cm'^, and T = 1.77 ±0.02, 
which are consistent with the results obtained using the XIS 
alone. 

The Suzaku XIS spectrum of 2RXP J 1301 at each epoch 
was then fit by using the cutoff power-law model with the 
parameters obtained above (namely, = 4.4 keV and ef = 
14 keV), to estimate the flux contamination of 2RXP J1301 
in each PIN spectrum. Spectral fitting was performed in an 
energy range of 2-10 keV after excluding an iron-K band 
of 6.2-6.8 keV in some occasions where we found possi- 
ble iron lines. The template model with the fixed spectral 
shape gives an acceptable fit in all cases, indicating the sta- 
bility of the spectral shape of 2RXP J 1301 from one epoch 
to another. The flux level of 2RXP J 1301 was also found to 
be stable, with the unabsorbed 1-10 keV flux in a range of 
F = (2.6-3.8) X 10"'' erg cm'^ s"'. 

We also fit Suzaku XIS spectrum of 2RXP J 1301 at each 
epoch by using an absorbed power-law model, instead of the 
cutoff power-law model of equation ([U, to quantify spectral 
stability. The single power-law model gives an acceptable fit 
in the 1-10 keV band in all cases, as summarized in Table[3] 
The observed spectral shape in this band looks similar to each 
other. The best-fit photon indices in epochs Szl-8 are found 
to be r ~ 1.2 ± 0.1, with the absorption column density of 
A^H = (2.7-3.0) X 10^^ cm"^. The spectral similarity from one 
epoch to another supports our assumption that the parameter 
sets of ec and ef obtained in Sz5h-6 are applicable to other 
epochs as well. 

3.4. Broadband Spectra 

Figures |6] shows the Suzaku XISh-PIN spectra of the 
PSR B 1259-63 system. Given the spectral similarity, two 
successive epochs were merged into a single observation, la- 
beled as SzlH-2, Sz3h-4, and Sz7h-8. For the PIN spectrum, the 
NXB and CXB as well as the conti-ibution from 2RXP J 1301 
were subtracted (see 33.3l l. We first performed spectral fitting 
using an absorbed single power law. The results are shown 
in Figure |6] and the best-fit parameters are summarized in Ta- 
ble g] 

The Sz3h-4 spectrum could not be well described by a sin- 
gle power law; there found wavy residuals in the XIS band, 
and, moreover, there exist flux deficits compared to the power- 
law model in the PIN spectra at high energies, indicating 
the presence of spectral steepening. It is interesting to note 
that similar wavy residuals to the best-fit power-law model 



TABLE 3 

Results of 2RXP J 1 30 1 Suzaku XIS Spectral Fitting 



ID A'h 


r 


Fi-io 


Xiid.o.t) 


10^2 cm-2 




lO^'^erg cm"2 §-1 





Szl 


2.78!<''| 


1 20+*-"' 


32.0!|;' 


1.09(192) 


Sz2 


2 77+0:25 


1.18±0.09 




0.91 (250) 


Sz3 


2 69-™ 22 


1.19±0.08 


36.4 ± 1.0 


0.95 (315) 


Sz4 


T qi+0.23 
^•^^-0.21 


1.20 ±0.08 




1.10(334) 


Sz5 


9 77+0.23 
^•"-0.22 


1 10+0.09 


-IT 7+LO 


1.02 (267) 


Sz6 




1.28 ±0.09 


26 4+" ' 


1.03 (267) 


Sz7 


1 17+0.09 

^•"-0.08 


31 4+'-'' 

-0.9 


1.05 (270) 


Sz8 


3 03+^'^^ 
-'■"-'-0.27 


1.26±0.10 


30.7+;;2 


1.00 (225) 



Note. — Fitting Suzaku XIS spectrum of 2RXP J1301 by a 
power law with photoelectric absorption in 1-10 keV. Absorb- 
ing column density A'h, photon index F, and unabsorbed 1-10 
keV flux Fi-10 are shown with 90% errors. 



are found in the XIS spectrum of the supernova remnant 
RX J 1 7 1 3 .7-3946, where the X-ray spectrum is indeed steep- 
ened ( Takahashi et al.l l2008l) . Imperfect subtraction of the 
2RXP J 1301 contamination is difficult to account for the flux 
deficit, because 2RXP J 1301 has a much steeper spectrum in 
the PIN band and hardly affects the PIN spectrum above 30 
keV. In fact, even if we artificially lower ef to 10 keV to in- 
crease the 2RXP J 1301 contamination, the presence of spec- 
tral steepening is still required. We conclude that the Sz3h-4 
spectrum exhibits steepening toward high energies. Also, the 
Sz1h-2 spectrum would require less significant steepening, the 
presence of which is not conclusive. The broadband Sz7h-8 
spectrum was found to be consistent with a simple power law 
with r= 1.63 ±0.02. 

We have proceeded to fit the Sz3h-4 spectrum with a bro- 
ken power-law model attenuated by photoelectric absorption, 
which yielded an acceptable fit (see Table |4|. The best-fit 
photon indices are Ti = 1.25!j5;[i^ ^2 = l-66!jj{]^, with a 
break energy of Sbr = 4.5!g 2 ^^^^ Sz3h-4 spec- 

trum with a broken power-law model restricting the energy 
band to 0.6-10 keV. In such an XlS-only case, we obtained 
the foflowing values: Ti = 1 .25 ± 0.04, T2 = 1 .64 ± 0.04, and 
Ebi = 4.4 ± 0.2 keV. Though the systematic error in the PIN 
spectrum itself would be large, its influence on the spectral 
parameters is therefore found to be small. Figure [T] presents 
the unfolded X-ray spectrum with the best-fit broken power- 
law model. For reference, the broken power-law model was 
applied also to the Sz1h-2 and Sz7h-8 spectra. As expected 
from the single power-law fitting, two photon indices for the 
Sz7h-8 spectrum are similar. 



3.5. Flux Upper Limit at 100 keV 

The HXD-GSO scintillation detectors cover the energy 
range of 50-600 keV. Since the contribution of 2RXP J 1301 
is expected to be negligible above 50 keV, the measurement 
by the GSO is important for testing the presence of spectral 
steepening inferred by the XIS -PIN fitting. There found no 
significant detection in any epoch with the HXD-GSO. We 
then set the upper limits on the source count rate as 2% of the 
background count rate, which corresponds to a ~ 3(7 limit. It 
should be noted that the uncertainty in GSO measurements 
is dominated by the systematic error associated with back- 
ground modeling. The upper limit of 1 x 10"'" erg cm"^ s"' 
was placed at 100 keV (in an energy range of 78-114 keV). 
The limit depends only weakly on the assumed photon index 
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TABLE 4 

Results of PSR B 1259-63 XIS+PIN Spectral Fitting 



ID 


Model A'h 


r, 






x'Cd.o.f.) 




10^2 




keV 







Szl+2 powerlaw 0.52±0.02 1.64±0.02 ■■■ ■■■ 1.16(604) 

Szl+2 broken PL 0.44!<55]2 L48!};5]^ 3.7±0.5 1.76±0.04 1.05(602) 

Sz3+4 powerlaw 0.53 ±6.02 1.43 ±0.02 ■■■ ■■■ 1.27 (500) 

Sz3+4 broken PL 0.43 ±0.03 l-25!}; j]2 4 5+0.5 l.66^°l 0.99(498) 

Sz7+8 powerlaw 0.48±0.02 1.63±0.02 ■■■ 1.12(490) 

Sz7+8 broken PL 0.47 ±0.02 1.59!}i ii^ 4 (fixed) 1.66!^ {i^ 1.11 (488) 

Note. — Fitting Suzaku XIS+PIN spectrum by a power-law or a broken power- 
law (PL) model with photoelectric absorption in 0.6-50 keV. Absorbing column 
density A'h, photon index Fi (and T2 in the case of broken PL), and the break 
energy Eb,- of a broken power-law, are shown with 90% errors. 



(F = 1 .5) for a reasonable range of F = 1-3. Since the averaged 
background level was almost unchanged during the monitor- 
ing campaign, differences in the upper limits from one epoch 
to another are negligible (only ^ 10%). Though the fact that 
the uncertainty comes from the systematic error prevents us 
from making a definitive statement about the significance, the 
upper limit set by the GSO depicted in Fig. |7] strengthens the 
presence of spectral steepening. 

4. INTERPRETATION OF THE SPECTRAL BREAK 

The most important finding with Suzaku would be the pres- 
ence of a remarkable spectral break by AF ~ 0.4 around 
£bi ^ 5 keV in the Sz3h-Sz4 spectrum. The low-energy part 
below £br can be described by a rather flat index of Fi = 
1.25 ±0.04. Here we consider the physical implication of 
this spectral structure within the framework of the shocked 
relativistic pulsar wind model, in which high-energy radia- 
tion comes from shock-accelerated electro ns and positrons 
pro duced at wind termination shock (e.g., iTavani & AronsI 
II997I) . Indeed, the synchrotron and IC radiation by the shock- 
accelerated pairs from the PSR B 1259-63 system as re- 
sults of the interaction between the relativistic pulsar wind and 
the Be star outfl ows have been studied extensively in previous 
work (Tavani & AronsI [l997t iKirk et alj [19991 iMurata et alj 
120031: tKhangulyan et_alj|2007!), on which our considerations 
are based. However, if the particles in the relativistic wind 
can be efficiently mixed with the surrounding dense medium, 
other radia tion channels such as hadronic interactions ma y be 
important dKawachi et al.l2004tlChernvakova et al.l2006l) : we 
do not consider such possibilities in this paper. Also, compre- 
hensive modeling of the overall phenomenology is deferred to 
a future pubUcation. 

4. 1 . Model of Compactified Pulsar Wind Nebula 

Here we describe our shock-powered emission model, 
which is broadly similar to what has been devised in previ- 
ous theoretical work mentioned above. The PSR B1259-63 
system can be regarded as a scale-down version of PWNe. 
The relativistic wind of the young pulsar is presumed to be 
confined by the stellar outflows of the form of wind and disk, 
leading to the formation of termination shock inside the pul- 
sar wind. The location of the termination shock should be 
close to the pulsar if the ram pressure of the stellar outflows 
as seen by the pulsar is large enough. The distance of the 
shock from the pulsar, r^, is assumed to be much smaller than 
the st ar-pulsar separation d: <ti d (e.g., iKhangulvan et al.l 
l2007h . This largely simplifies the calculations below. For ex- 



ample, the stellar photon density at the wind shock region as 
well as the IC scattering angle for the observed gamma-rays 
can be approximated as constant at a given phase. 

The shocked relativistic wind, at the distan ce of from the 
pulsa r, is magnetized with strength of (e.g., iTavani & AronsI 
Il997h : 




where is the spindown power of the pulsar, and a is the 
magnetization factor defined by the ratio of the Poynting and 
kinetic energy flux in the pulsar wind. A factor of 3 in equa- 
tion ^ is introduced to account for the shock compression, 
and (7 <C 1 is assumed in obtaining equation ([3]). The spin- 
down power is measured to be /ip ~ 8 x 10^^ e rg s~' . The mag- 
netiza tion factor was adopted as cr = 0.02 in ITavan i & AronsI 
( 119971) . The Ken nel & Coroniti model of t he Crab Nebula 
yielded a = 0.003 dKennel & Coronitil l 1984a'), which was de- 
termined to satisfy the flow and pressure boundary conditions 
at the outer edge of the nebula. 

The radiation energy density at a distance d from the Be 
star can be estimated as: 

where is the luminosity of the Be star: ~ 7 x 10^' erg s"' 
for Teff = 27000 K and R^, = 6Rq. (L* is somewhat uncertain. 
For a larger stellar radius of = IQRq, becomes larger by 
a factor of 3.) The stellar light dominates over the magnetic 
field of the pulsar wind at the termination shock. However, 
because the IC scattering of TeV electrons on stellar photons 
occurs in the deep Klein-Nishina regime, the IC emission can 
be largely suppressed. The synchrotron cooling is expected 
to overtake the IC cooling at high energies. The larger lumi- 
nosity of the observed synchrotron X-ray emission compared 
with the IC gamma-rays during periastron passage suggests 
that the synchrotron cooling is indeed faster than the Klein- 
Nishina IC cooling at TeV energies. 

Given the postshock flow speed of ~ c/3 in the pul- 
sar wind, the adiabatic loss can compete with the radiative 
losses. The adiabatic loss rate is adopted as = -7/fad 
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FIG. 6.— Suzaku XIS+PIN spectrum of PSR B 1259-63 in the 0.6- 
50 keV band and the ratio of the data and the best-fit absorbed power- 
law model: (a) Szl-l-2, (b) Sz3-l-4, and (c) Sz7-l-8. The contaminating 
flux from 2RXP J 1301 in the PIN spectrum was subtracted in the 
way described in 33.31 Superposed on the spectral data are the ab- 
sorbed single power-law models, which do not provide a good fit to 
the Szl-l-2 and Sz3-l-4 spectra. 

with fad = C^s/c. The factor ^, likely ^ ~ 3-10, is treated 
as a free parameter In Figure [8] we show the energy- 
dependence of the cooling timescales both for radiative (syn- 
chrotron and IC) and non-radiative (adiabatic loss) losses, at 
r = 30 days, where the importance of adiabatic losses is il- 
lustrated. Recently, it has been shown that the hydrodynam- 
ics of the interaction of the pulsar and stellar winds do pre- 
dict the predom inance of the ad iabatic loss over the radiative 
losses ( Khangulvan et alj|2008t) . The discussion on the cool- 
ing regimes in the PSR B 1259-63 system can be found in 
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Fig. 7.— Broadband unfolded X-ray spectra of PSR B 1259-63 
obtained by using the Suzaku data in epochs Sz3 and Sz4. The low- 
energy absorption is corrected for by the best-fit model. The best-fit 
broken power-law function is overlaid. 

iTavani. Arons. & Kaspil ( 119941) . 

The Lorentz factor of the relativistic wind is denoted by 71 ; 
the unshocked pulsar wind is presumed to contain with 
an ener gy of ^inigC^. Using a = 0.003, iKennel & Coronitil 
(Il984bl) obtained 71 = 10^ for the Crab pulsar to account for 
the observed nonthermal radiation from the nebula. 

Unshocked cold pairs in the relativistic wind are as- 
sumed to be accelerated at the termination shock and be in- 
jected into the downstream postshock flow. The energy dis- 
tribution of the injected nonthermal pair is described by 
2(7) = 2o7~''exp(-7/7ni) for 7 > 71. A low-energy cutoff 
with a step function is assum ed: Qil) = for 7 < 71. Re- 
sults of simulation by Hoshino et al.l d 19921) for the relativistic 
coUisionless shock in electron-positron-ion plasma suggest a 
power-law spectrum with /7 ~ 2 can be formed above 71 , and 
a fraction e ~ 0.2 of the flow energy can be transfeiTed to the 
nonthermal pairs, under the condition that the upstream 
flow energy carried by ions exceeds that of pairs. We note 
that 7m should be limited by rapid cooling of the pairs at 
highest ene rgies. The OSSE r esults in the year 1994 indicated 
7m > IO71 (iGrove et al.ll995h . The normalization Qq is deter- 
mined by the spindown power of the pulsar and the parameter 
e, through the relation: 



eE„ 



(6) 



The energy distribution of the accelerated pairs inte- 
grated over the entire radiation zone, A^(7), can be derived 
by the well-known kinetic equation. We consider simplified 
cases where both the magnetic field and radiation field do not 
change within the radiating zone. The cooling timescale of 
the accelerated pairs is much shorter than the dynamical 
timescale of the binary. Therefore, the energy distribution 
N((j),j) at a given orbital phase (f) is reduc ed to be a steady- 
state solution (e.g.. lKhangulyan et alJl2007l) : 



1 



(7) 



where b((j>,j) = -(7syn + 7ic + 7ad) is the energy loss rate as 
given by a sum of the synchrotron, IC, and adiabatic loss rates. 
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Fig. 8. — Cooling timescales as a function of electron (positron) 
energy E = ^nieC^ at r = 30 days witii tlie following parameters: B = 
0.4 G, d = 1 .8 AU, hi = lOn/c, and n = O.lSd. 



TABLE 5 

Parameters of Synchrotron-IC models 









Fixed 






Free 


T 


e 


a 


71 




P 


Em 


c 
















(TeV) 






—15 days 


0.1 


0.01 


4x 


105 


1.9 


10 


0.05 


3 


+30 days ■ 


0.1 


0.01 


4x 


105 


1.9 


10 


0.15 


10 


+618 days 


0.1 


0.01 


4x 


105 


1.9 


10 


0.50 


2 



Note. — e (a fraction of the spin down power channelled into 
the accelerated pairs), a (the magnetization factor of the pulsar 
wind), 71 (the Lorentz factor of the pulsar wind), p (the acceler- 
ation index), E,„ = ^,^mgC^ (the maximum energy of accelerated 
pairs), C, = rs/d (the distance of the termination shock from the 
pulsar divided by the pulsar-Be star separation), and ^ (the param- 
eter to describe the adiabatic loss rate). 

We are now ready to calculate the synchrotron and IC 
spectra using Nijp.,^), and compare them with the obser- 
vations. Anisotropic IC scattering of stellar photons is 
taken into account by utilizing the IC kernel derived by 
|Ah aronian & Atovan ( 1981), as the IC gamma-ray emission 
depends largely on the a ngle between the star-pulsar an d 
pulsar-observer directions (iKhangulvan & Aharonianll2005h . 
We assumed a stellar radiation field to be radial rays originat- 
ing from the point source, the Be star The effect of finite size 
of the star does not affect the resultant radiation spectrum. 
However, the radiation field can be modified appreciably by 
the presence of the equatorial disk, which has to be taken into 
account by future work to model the TeV gamma-ray spec- 
tra. We calculate the attenuation of TeV gamma-rays due to 
pair production with stellar photons, though the 77 opacity 
is small, typically txp(-T^^) ^ 0.9, for our choice of the op- 
tical luminosity. Whereas more significant absorption, up to 
expC-T^^) ^ 0.6, can be expected if we adopt larger stellar 
luminosity (lDubusll2006l) , this does not affect our main con- 
clusions. 

4.2. Applications to the Observed Spectra 

To describe the broadband emission of the compactified 
PWN in variable environments, two parameters, C, = rjcl and 
^, are allowed to vary along the orbit. We note that C,, which 
is controlled by ram pressure of stellar wind/disk, determines 



magnetic field strength. The two parameters describe physi- 
cal conditions external to the pulsar wind. The fact that we 
observed a spectral break during the disk crossing suggests 
that the external conditions indeed regulate the position of 
the break. The other parameters are frozen to be the fol- 
lowing values: e = 0.1, a = 0.01, 71 = 4 x 10^, p = 1.9 and 
Em =10 TeV. They are broadly similar to what have been 
previously obtained in the literature mentioned above, though 
we made some fine-tuning; 71 is adjusted to explain the spec- 
tral break, and p is chosen to match with the observed X-ray 
spectral slope in Sz7h-8 (see below). The uncertainties of the 
fixed parameters are discussed in ^4.31 The physical parame- 
ters are summarized in Table |5] 

We first reproduce the Sz7h-8 spectrum, which is described 
by a simple power law of F = 1.63 ± 0.02 in a wide energy 
range of 0.6-50 keV, a long with the HESS spectrum in 2004 
jAharonian et al.ll2005l) . Since our treatment of the radiation 
field ignores the presence of the equatorial disk, we require 
only a crude match with the average flux level reported by 
HESS. (The TeV gamma-ray flux changed by a factor of 2- 
3 from one month to the next.) Nevertheless, the TeV flux 
and shape put meaningful constraints on the model. Using 
the configuration of the binary system at r = +30 days, we 
found that the following parameters give a reasonable fit to 
the X-ray and gamma-ray spectra: C = 0.15 (r, = 0.27 AU and 
B = 0.38 G), and ^ = 10. In Figure|9] the spectral energy dis- 
tribution (SED) of the synchrotron and IC radiation is drawn. 
The cooling of the accelerated responsible for the soft X- 
ray emission is dominated by adiabatic loss, which ensures 
the hard X-ray spectrum. The shape of the TeV gamma-ray 
spectrum obtained with HESS can be reproduced by this mod- 
eling as well. 

We then proceed to model the Sz3+4 spectrum, where a 
spectral break was observed. We relate the X-ray spectral 
break to the low-energy cutoff in the electron injection func- 
tion, which can be ascribed to the Lorentz factor of the pul- 
sar wind 7i. In Figure |9] we present calculation of the 
synchrotron and IC models at r = -15 days with ( = 0.05 
(r, = 0.057 AU and B = 1.8 G), and ^ = 3. The X-ray spec- 
tral shape, the very flat spectrum below Cbr in particular, is 
nicely reproduced by this model. With this scenario, the ap- 
pearance of the break in the X-ray band is due to the enhanced 
magnetic field. Our model predicts lower flux in the TeV band 
at the time of a higher break energy. 

The observed break energy is related to the Lorentz factor 
of the pulsar wind: 

[3 eh , 

^b'=\/o 2^7? (8) 

^4('-^V-lL_VkeV. (9) 
V1.8Gy V4xlOV ^ ^ 

It is interesting to note that the X-ray spectrum below ebi- is 
predicted to have a universal shape as long as the cooling is 
dominated by adiabatic loss and the energy distribution 2(7) 
below the low-energy cutoff is harder than 7"^ Given the 
adiabatic-loss dominance (b cx 7), equation O yields 

A^(7) cx 7"' for 7 < 71 , (10) 

irrespective of acceleration index p. This form of the energy 
distribution of the cooled pair population is translated into 
a synchrotron power-law spectrum with photon index F = 1 . 
Therefore, the synchrotron spectrum has a universal form with 
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Fig. 9. — Broadband spectra of the PSR B 1259-63 system along with the synchrotron (X-ray) and IC (gamma-ray) models. The average 
HESS spectrum in 2004 is indicated to allow a comparison with the IC calculation. The flux upper limit at 200 MeV is obtained with the 
EGRET observations in 1994. Red: Synchrotron-IC model at r = 30 days in comparison with Suzaku X-ray spectral data for Sz7-l-8. Blue: 
r = -15 days with Suz aku Sz3-(-4 which exhibits a spectral break. Green: r = 618 days (apastron) with AxSCA X-ray spectrum around apastron 
( lHiravamaetaUl999h . 



r = 1 for e <C Ebr- In fact, the result of the Suzaku observa- 
tions, namely F ~ 1 .2 immediately below ebi , agrees well with 
what theoretically expected. 

Let us check if this model can accommodate the X-ray 
observations around apastron, where the physical conditions 
change dramatically from periastron. We found the follow- 
ing parameters are compatible with the X-ray observations at 
apastron: ( = 0.5 (r^ = 4.7 AU and B = 0.02 G), and ^ = 2. 
(A condition of ^ <C 1 is not well satisfied, so our calculations 
need a revision if we aim at performing precise modeling. ) In 
Figure |9] we show the ASCA spectrum and the model appro- 
priate for apastron, r = 618 day. The effective photon index of 
the model is Fsyn = 1 .7, which is in agreemen t with the X-ray 
spectrum, F = 1 .6 ± 0.2, obtained with ASCA (iHirayama et alJ 
[1999). 

4.3. Basic Properties of the Pulsar Wind of PSR B1259—63 

One of the most fundamental parameters of a pulsar wind, 
the Lorentz factor of the wind 71, can be estimated by our 
modeling in the previous section through the identification of 
the spectral break with the low-energy cutoff of the acceler- 
ated pairs. We adopted 71 = 4 x 10^ to account for the ob- 
served break position. It depends only on the square root of 
postshock magnetic field, cx \fB, which can be constrained 
well by the X-ray to TeV gamma-ray flux ratio. Therefore, 71 
is hard to be altered more than by a factor of few. 

It should be noted that the wind Lorentz factor will pos- 
sibly be measured by the Fermi Gamma-ray Space Tele- 
scope during the next periastron passage. The Comptoniza- 
tion of the unshocked pulsar wind inevitably leads to the 
formation of an additional bump-like gamma-ray spectrum 
at £ 7imeC^ (assuming the Klein-Nishina limit), which 
overwhelms the IC component in Fig. |9]in the 0.1-10 GeV 
band ( Khangulvan et al.ll2007l) . The X-ray and GeV-to-TeV 
gamma-ray data during the next periastron passage will be 
crucial to test the models of the PSR B 12 59-63 system, and 
the lon g-standing pulsar wind paradigm jKennel & Coronitil 
Il984ah in general. 



The other fundamental parameter in the standard pulsar 
wind paradigm is the magnetization factor a. While the 
magnetic field strength can be well constrained by the X-ray 
and TeV flux ratio as B ~ 0.5 G during the disk crossings, 
the (T-parameter itself is allowed to vary in a wide range of 
10"'* < (T < 1 depending on r,. The condition r, < d required 
to explain the spectral variability coiTesponds to < 1, and 
the energetics requirement of e < 1 (together with ^ ^ 10) 
gives a lower bound as cr > 10"'*. 

The magnetic field of Z? ~ 0.5 G constrains particle acceler- 
ation mechanism(s) at work in the termination shock region. 
Synchrotron cooling time of 10 TeV electrons is about 100 
sec. To overcome the fast cooling, the acceleration timescale 
of face ~ 100 sec ^ 50r]^/c is necessary at £ = 10 TeV, where 
rr, de notes the Larmor radius of pairs (see iKhanguIyan et alJ 
2007). According to the magnetosonic acceleration model of 
iHoshino et al.l (Il992h where the magnetosonic waves collec- 
tively emitted by the shock-reflected protons (or heavy ions) 
are resonantly absorbed by pairs, the inverse of the gyrofre- 
quency of protons with a Lorentz factor of 71 determines the 
acceleration timescale; interestingly, face ^ '«pC7i/(eB) ^ 50 
sec is compatible with the requirement posed by our model- 
ing. 

The efficiency with which the pulsar spindown power is 
transferred into the nonthermal e^pairs, namely e that is de- 
fined by equation (|6]), was adopted as e = 0.1 in our mod- 
eling, which i s also consistent with the numerical result of 
iHoshino et al.l (Il992l) . It should be noted, however, that the 
determination of e would be subject to the degree of relativis- 
tic Doppler boosting if possible rela tivistic motion of the post- 
shock flow (Bogoval ov et alj|2008l) plays a major role in the 
nonthermal emission. 

4.4. Alternative Model for the Spectral Break 

We have explored another possibility to explain the X-ray 
spectral break — the observed break may be accounted for by 
the transition of main cooling mechanisms; below the break, 
cooling of pairs is assumed to be IC cooling in the Klein- 
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Nishina regime, while synchrotron cooling takes over above 
the break. Due to the inefficient IC cooling in the Klein- 
Nishina, the electron energy distribution can be much harder 
in a certain energy band, as compared to the energy bands 
wherein the IC cooling in the Thomson regime or synchrotron 
cooling are dominant. With this interpretation, ^ > 100 is nec- 
essary to make the adiabatic loss unimportant. 

To match the flat X-ray spectrum below Sbi, the accelera- 
tion index should be very flat, p ~ 1.0. However, the very 
flat spectrum of accelerated would be inconsistent with 
the conventional particle acceleration model as well as the 
observations of PWNe on parsec scales. Also, the signifi- 
cant reduction of the X-ray flux around apastron requires a 
drastic change of ^ (say, from f = 100 to ^ = 3). In view of 
these uncomfortable requirements, we prefer the adiabatic- 
loss-dominated model in ^4.21 where the X-ray break is as- 
cribed to the appearance of the Lorentz factor of the relativis- 
tic pulsar wind. 

5. CONCLUSIONS 

We have presented results of new X-ray observations of the 
pulsar-Be star binary system PSR B 1259-63 performed eight 
times near periastron using the Suzaku satellite. The moni- 
toring observations with Suzaku, coupled with previous data, 
allow us to identify the two bump in the X-ray lightcurve 
around periastron, which can be ascribed to the pulsar's en- 
trance to the equatorial disk of the Be star. In most cases, 
the X-ray spectra can be well fit by a simple power law, as 
observed in previous observations. However, we found evi- 
dence of a spectral break when the X-ray spectral shape below 
10 keV is described formally by a rather hard power law with 
photon index T ^ 1.3. Applying a broken power-law model 
to the Sz3h-4 spectrum in the 0.6-50 keV band, we found 
that the index changes by AT ~ 0.4, from Ti = 1.25^o;o| to 
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= 1 .66t!^ °ol, around et,- = 4.5!!]:^ keV. 
We considered a compactified PWN model, in which the X- 
ray emission is accounted for by shock-accelerated relativis- 
tic pairs via synchrotron process. The spectral properties 
of the X-ray and TeV gamma-ray emission are found to be 
reconcilable in a compactified PWN model whose parameters 
are adjusted within a reasonable range. We ascribed the spec- 
tral break to the low-energy cutoff of the accelerated pairs, 
which is supposed to be linked directly to the wind Lorentz 
factor We here emphasize that our choice of parameters (see 
Table |5]l is definitely not the unique one that can fit the ob- 
served spectra. Nevertheless, the value of the wind Lorentz 
factor, 7i ~ 4 x 10^, determined by the observed break energy 
in the X-ray band is robust within a factor of few, because it 
depends only on the square root of postshock magnetic field, 
which in turn can be deduced from the X-ray to TeV gamma- 
ray flux ratio as B = 0.4-2 G during the periastron passage. 
Also, the magnetic field estimate indicates the magnetization 
factor should be in a range of 10"'* < cr < 1. The determi- 
nation of the wind Lorentz factor through the X-ray spectral 
break, if confirmed, provides a direct probe of the standard 
pulsar wind paradigm in which pulsar's spindown power is 
supposed to be carried off by a relativistic wind with Lorentz 
factor 7i > 10"^. 
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